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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

BLOCKAGE CORRECTIONS FOR THREE-DIMENSIONAL-FLOW CLOSED- 
THROAT WIND' TUNNELS, WITH CONSIDERATION 
OF THE EFFECT OF COMPRESSIBILITY 
By John G. Harriot 


SUMMARY 

Theoretical blockage corrections are presented for a body of 
revolution and for a three-dimensional unswept wing in a circular or 
rectangular wind tunnel. The theory takes account of the effects 
of the. wake and of the compressibility of the fluid, and is based 
on the assumption that the dimensions of the model are small in 
comparison with those of the, tunnel throat. Formulas are given for 
correcting a number of the quantities, such as dynamic pressure and 
Mach number, measured in wind— tunnel tests. The report presents' a 
summary and unification of the existing literature on the subject." 


INTRODUCTION 

When a model is placed in a closed— throat wind tunnel there is 
an effective constriction or blockage of the flow at the throat of 
the tunnel. The effect of this blockage is to increase the velocity 
of the fluid flowing past the model; if the model is not too large 
relative to the tunnel throat, this veloqity increment is approxi- 
mately the same at all points of the model so that the model is 
effectively working in a uniform stream of, fluid the velocity of 
which is, however, greater than the free-stheam velocity observed 
at some distance upstream of the model. It is therefore necessary 
to correct the observed velocity, dynamic pressure, Mach number, 
and other measured quantities for the effect of this constriction. 
This correction is frequently called the correction for "solid 
blockage.’ In addition to this solid— blockage correction, a cor- 
rection for "wake blockage" is also necessary if the true dyn ami c 
pressure and Mach number at the model are to be determined." This 
wake blockage arises because the fluid is slowed down in the wake 
and consequently must be speeded up outside the wake. A further 
effect of the wake is to produce a pressure gradient which must be 



2 


MCA RM.No. A7B28 


considered in correcting the, drag coefficient. 

Formulas for the solid-blockage correction for a model mounted 
in a two-dimensional— flow wind tunnel are given in references 1, 2, 
3, and 4. The first order effects of the compressibility of the 
fluid on these corrections are given in reference as well as in 
references 3 and k. References "2, 3 , and 4 consider also the wake— 

correction, and reference 3 also considers the drag cor 

rection due to the pressure gradient caused by the wake. The 

formulas given in reference 3 for the solid— and wake— blockage cor 

recti ons will usually be found most convenient whenever the engineer 
is confronted by a practical problem of determining the corrections 
for any configuration met in his experimental work. 

The solid-blockage correction for a model mounted in a three- 
dimensional— flow wind tunnel has, been given in a number of different 
forms by different authors. Not only do different authors give the 
correction for the same configuration in different forms but no one 
author gives formulas which are applicable to both fuselages and 
wings in tunnels of various shapes; for thisreason, the engineer 
confronted with a correction problem may have to refer to several 
reports to get the complete solution of his problem. Moreover the 
modifications of the formulas for the first order effects of fluid 
compressibility are given incorrectly in some cases. References 
1 and 2 give a formula for the solid— blockage correction for a body 
of revolution in a circular or rectangular tunnel for the case of 
incompressible flow. In references 5 and 6 the effect of the com- 
pressibility of the fluid on this correction is discussed, but the 
result given is incorrect.' Reference k gives a formula for the 
solid-blockage correction for a body of revolution and for a three- 
dimensional wing in the 7- by 10-foot wind tunnel. The modification 
for compressibility is correct for the wing but' wrong for the body 
of revolution. Reference 7 gives a formula for the solid-blockage 
correction for a body of revolution in a circular tunnel, correctly 
taking account of the effect of the compressibility of the fluid. 
References 8, 9, and 10 give a formula for the solid-blockage cor- 
rection for any body in a circular wind tunnel together with the 
appropriate constants for a body of revolution and for a rectangu- 
lar wing having various spanr-to-diameter ratios; the modification 
of this formula to take account of compressibility is correctly 
given. 

It is clear that no one report gives all the necessary for- 
mulas together with the appropriate constants and compressibility 
modifications to enable the engineer to calculate the solid- 
blockage correction for any case with which he may be confronted. 
Moreover, when the re suits:: of two or more reports overlap, the forms 
are frequently different so that it is not obvious whether the 
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results are In agreement. It is the purpose of the present report 
to summarize and extend the .results of the previously mentioned 
reports. Formulas are given for the' calculation of the solid— 
"blockage correction for a "body of revolution or a three-dimensional 
unswept wing in a circular or rectangular tunnel. These formulas 
contain two constants, one depending on. the shape of the body and 
the other on the shape of the tunnel and the ratio of wing span to 
tunnel breadth. (This ratio may be taken to be zero for a body of 
revolution.) Values of the first constant for various bodies of 
revolution and for a number of frequently encountered wing— profile 
sections are given. Values of the second constant for a circular 
tunnel and for rectangular tunnels of a number of commonly encoun- 
tered breadth—' to-height ratios are given for various wing— span-to— 
tunnel— breadth ratios. Some of these values have been taken from 
the previously mentioned reports; whereas others appear for the 
first time in the present report. The discussion is limited to 
bodies centrally located in the wind tunnel. 

The wake— blockage correction for a model mounted in a three- 
dimensional— flow tunnel is given in references 4, 8, 9, and 10. 

For the case of incompressible flow the formulas are in agreement, 
but the modification to .take account of compressibility given in 
reference 4 differs from that given in references 8, 9 , and 10. 

This matter is discussed in the present report. In reference 'll 
the corrections for the pressure gradient due to the wake are 
given with the correct compressibility factors. 

All the final correction formulas together, with directions for 
their use are given in the final section entitled "Concluding 
Bemarks." Mathematical symbols are defined as introduced in the 
text. For reference, a list of the more important symbols and their 
definitions is given in appendix B. 


SOLID BLOCKAGE IN INCOMPRESSIBLE FLOW 

In studying the flow over a thin airfoil of small camber at a 
small angle of attack it has been shown that the effects of camber 
and thickness may be considered independently. In treating the 
problem of wall interference, it is again convenient to consider 
the thickness and camber .effects separately. The camber effect, as 
pointed .out in reference 3., contributes nothing to the blockage 
correction. Consequently it suffices to determine the blockage 
correction for symmetrical bodies at zero angle of attack. This 
means that for wings it is necessary to consider only the base pro- 
file of the airfoil, the base profile being defined as the profile 
the airfoil would have if the camber were removed and the resulting 
symmetrical airfoil placed at zero angle of attack; bodies of 
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revolution need be considered only at zero angle of attack. 


.Rectangular Tunnel 

Three-dimensional wing .— The blockage correction for a three- 
dimensional wing in a rectangular tunnel is considered in refer— - 
ence 4. Numerical values are given only for a wing of 6-foot span 
in a 7 — t»y 10-foot wind tunnel. The method may, however, be applied 
to any rectangular tunnel and any span-to— breadth ratio. 

Consider a rectangular tunnel of height ' H and breadth B. 

Suppose that the wing span 2s is. in the B-direction so that 2s/B 
is the span-to-breadth ratio. As in reference 4, let the wing be 
represented by a series of finite lines of sources and sinks. (See 
fig. 1. ) The strengths of these sources and sinks are assumed to 
depend on the airfoil profile and the determination of these strengths, 
which is a two-dimensional problem, is explained presently. If each 
line retains. the same strength from one end to the other, the wing 
section cannot be exactly constant. It will thin down at the 
extreme tip and the plan form will not be exactly rectangular, but 
neither of these features is such as .to detract from its usefulness 
for the present purpose > which is to represent an a<3tual wing 
sufficiently well to enable a calculation to be made of the velocity 
along the tunnel axis "induced" by images of the wing. It is this 
velocity induced by the images which represents the effect of the 
tunnel walls on the, velocity at the model, and which is the solid- 
blockage correction. 


Consider the image line source CD of strength Q per unit 
length. The veloci.ty potential of a three-dimensional source of 
strength Q&y (volume per unit time) is — Q,6y/4jtr. It follows that 
the component velocity along the tunnel center line induced at A 
by the source element Q5y is 


8UA “ ^ &y 


where g and 
r 2 « (mB-y) 2 + n 2 ^ 


are defined in figure 1. Putting 
+ g and integrating from — s to 


s gives. 


AUa Wn^+g^) [ /“ 


mB + 8 


mB — s 


/ n i2 H 2 +g 2 +(mB+s.) 2 v/ n 2 H 2 +g 2 + (mB-s ) 2 _■ 
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for the single line source CD. Now CD is one image of one finite 
line source used to represent the wing in the tunnel. In order to 
find the velocity induced "by all the. images of this particular line 
source, it is necessary to add the results obtained from equation (l) 
by giving m and n all positive and negative integral values 
except m =* n = 0. Thus for a single line source. 


U A 


=: 9E 
4rr 


r 


n 2 H 2 +g 2 


niB 


+ s 


mB — s 


L ^n 2 H 2 -fg 2 * ( mB + s ) 2 A i 2 H 2 +^+( mB-s ) 2 J 


( 2 ) 


where the prime denotes that the term m = n = 0 is to be omitted 
from the summation. Equation .(£?) may be rewritten 


u A = a r E s B 

A 2jtH 3 Ih B H 


where 


a 


£ s B\ = B_ 
H B* H J 2b 


1 


1 _ 
n 2 +(g/H) 2 


m + s/B 

- / n^+(g/H) 2 +(m+s/B) 2 (B/H) 2 


m - b/E 

^ n 2 + ( g/H ) 2 + (mr-s/B ) 2 (B/H) 2 - 


( 3 ) 


It is convenient to define 


so that 


{ g 

s 

B N 

1/ 


E 3 

Hy 

' 2 \ 


_B\3/2 

nHy 


o 


f 8 8 B \ 
\E 3 3’ eJ 


UA = — Ss§. - *1/2 
A (EE) 3 / 2 



s b\ 
B* H J 


(b) 


Now a and r depend only very slightly on g/H and so it usually 
suffices to take g/H =0 in the evaluation of these quantities. 

Any wing profile can be represented by a suitable distribution of 
sources and sinks along the chord. It follows that the total 

induced velocity due to the wing images is obtained by summing over 
this distribution and is approximately ' 
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(5) 


It may "be noted here that setting g/H = 0 in the evaluation 
of a and t is equivalent to representing the wing by a line 
doublet of strength EQg (analogous to references 1 and 2) instead 
of by a distribution of sourees and sinks. 

The quantity EQg of equation (5) can be determined approxi- 
mately from the wing profile and this determination is a two- 
dimensional problem. From reference 2, but with the notation' of 
reference 3> there is obtained 


EQg = | Ac 2 U' 


( 6 ) 


where 

c airfoil chord 

U’ apparent free— stream velocity at airfoil as determined from 
measurements taken at- a point far ahead of model 

A a factor dependent on shape of base profile • 

Substitution of equation (6) into equation (5) yields 

AlU’ _ 2sc a £ 3/2 At f q s B \ 

U* '(BH)3/2 16" \ ’ B’ K J 


2set n 3 / 2 , A -f (\ • s. B 'X 
C 3 /2 16 t/c \ 3 B 3 H ) 


(T) 


where 

t maximum thickness of airfoil. ' 

C cross-sectional area of tunnel 

\ ' 

If V denotes the volume of the wing, then V = 2sctki where 
Ki 5 1 depends on the shape of the base profile. Equation (7) may 
be rewritten 
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or 


where 


Aj.tr* _ KifV 
U» " C 3 / 2 


( 8 ) 


AiU* _ K2 T 2sct 

~ C 572 — ' 


( 9 ) 




( 10 ) 


K x 


* 3/Z A 1_ 
16 t/c 


( 11 ) 


, Ks 


TT 3 / 2 A 

16 t/c 


( 12 ) 


Xt is clear that ' t depends only on the tunnel shape and the 
wing— span— to-tunnel— breadth ratio; whereas K lt K 2 depend only on 
the shape of the base profile. 


Ihe factor A can be determined for any base profile from the 
relation (references 2 and 3) 


where 


A = ~ J' 1— P V l+(dy^./dx) 2 d^2 'j 


( 13 ) 


y t ordinate of base profile at chordwise station x>: 

<*y t /cbc slope of surface of base profile at x 

P base— profile pressure coefficient at x in an 

incompressible flow 


Values of A for a number of base profiles are given in reference 3. 
Thus the value of K 2 can be calculated from equation (12), The value 
of Kj. is immediately found from the area of the base profile which 
may be calculated, for example, by a numerical integration from the 
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ordinates of the base profile. As soon as' ki is known, Ki can 
then be calculated from equation (ll) . The evaluation of t 
requires the summation of the infinite series in equation ( 3 ). The 
summation of this series is explained in appendix A. 

Values of r for rectangular tunnels of various breadth— to- 
height ratios and for various wing— span— to— tunnel— breadth ratios 
are given in table I and figure 2. Values of Ki and K 2 for 
various base profiles are given in tables II and III. . The 'choice 
between the two formulas (8)- end ( 9 ) is entirely a matter of con- 
venience and should be decided in the light of the available data. 

Body of revolution.— The blockage correction for a body of 
revolution in a rectangular tunnel is considered in references 1, 2, 
and 4. In reference 4 numerical values are given for some average 
streamline body of revolution in a 7— by 10-foot wind tunnel; 
whereas in references 1 and 2 numerical values are given for prolate 
spheroids and Rankine Ovoids in square and duplex (breadth equal to 
twice the height) wind tunnels. Either the method of reference 4 
in which the body of revolution is represented by a suitable dis- 
tribution of sources and sinks along its chord or the method of 
references 1 and 2 in which the body is represented by a. doublet of 
suitable strength at its center may be used tp obtain results for 
any streamline body of revolution in any rectangular tunnel. Since 
the method of reference 4 was used for the three-dimensional wing, 
it is instructive to use the doublet method of references 1 and 2 
for the body of revolution case, although both methods give the 
same results. 

Consider a body of revolution of maximum thickness t and 
length c centrally located in a rectangular tunnel of breadth B 
and height H. As in reference 2 the body may be represented by a 
doublet of strength p given by the equation 


p = £XtV (14) 

4 

where X. is a constant depending only on the shape ana fineness 
ratio of the body. The velocity induced at the model by the tunnel 
walls is the same as that induced by a doubly infinite array of 
images of the doublet and is given by 


(n^+m^B 2 ) 3/a 


AiU 1 


t 


(15) 
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the summation being taken over all positive and negative integral 
values of m and n except m = n * 0. If g/H = 0 equation (3) 
may be rewritten 


( 0 , s , s 

)= ~ ^ 

— 1 * 

J, 

' (m+s/BjF^ - (m-s/BjE^' 

\ B Es 

' 2s L 

n 2 

®mnFmn 



. 2m 

EmnFmn CU+a/BjF^+fm-s/BjEjnJ, 


(16) 


where the quantities Ejjj and F^, which are Introduced for con- 
venience, are defined by the equations 


E mn = ' / n 2 +(m+s/B) 2 ‘(B/H) 2 


Fmn v/ n 2 +(m-s/B) 2 (B/H) 2 


It follows that 



1 

[n 2 + m 2 (B/H) 2 ]3/2 


(17) 


Substitution of equation (17) into equation (15) yields 


AxU* = — H 


hitH 3 


a yO, 


0 , z) 

H J 


If 4 is replaced by its value from equation ( 14) and equation (4) 
is used, there is obtained 


A X U’ 

U* 


Tt3//2 *-t 3 r /' q q B ct 2 it 3 / 2 Xt 

8 (BH) 3/2 ^ c 3/2 8 C 


t(o, 0, - 
\ H 


(IS) 
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If V denotes the volume of the body of revolution* then 
V ±= K 2 ct 2 where K 2 — ir/4 depends on the shape of the meridian 
section of the body. (k 2 = tt /4 for a right circular cylinder 
whose meridian section is. clearly a rectangle.) Equation (l 8 ) nay 
be rewritten 


A 3 U 1 _ K 3 tV 

U' C 3 / 2 . 


or 

AxU 1 _ g 4 rct 2 
U» “ C 3 /s 

where 

t > t (°, 0 , |) 

1 K = 1_ 

3 . 8 C K 2 

„ It 3 /2 Xt 


(19) 


( 20 ) 


( 21 ) 

( 22 ) 

(23) 


It should be noted that r for the case of the body of revo- 
lution is the same as for the limiting case of a wing yhe n the span 
approaches zero. As pointed- out in reference 2, X may be calculated 
for any body of revolution whose pressure distribution is known. The 
necessary formula is 

x k (f) J (j) ^1-P /i+(dy/dx) 2 d(|) (24) 


where 

y radius of body at chordwise station 

dy/dx slop© of meridian section at x 


x 
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P pressure coefficient at x in incompressible flow 

Values of X for prolate spheroids and Bankine Ovoids are given in 
references 1 and 2. As soon as X is known for a body, K 4 can 
be calculated at once from equation (23).. The value of k g may' be found 
from the volume of the body of revolution which may be calculated 
for example by a numerical integration, and K 3 can. be -calculated from 
equation (22) as soon as ks is known. 

Values of j for rectangular tunnels of various breadth— to— 
height ratios are given in table I and figure 2. Values of K 3 
and K4 for various bodies of revolution are given in tables IV 
and V. Some of these bodies are drawn in figure 3 • Unfortunately 
these tables are rather incomplete and moreover fuselage shapes 
used in practice are extremely varied. However, in table IV it is 
observed that the values of K 3 do not depend very strongly on the 
shape of the body, although they do depend on the thickness ratio. 

For this reason it appears that for most fuselages it will be 
sufficiently accurate to use the values of E 3 given for the 
NACA 111 bodies. As the values of K 4 given in table V are more 
dependent on the body shape it is recommended that equation (19) 
and table IV be used in preference to equation (20) and table V 
whenever possible. 


Circular Tunnel 

Body of revolution — It is convenient to consider the case of a 
body of revolution before considering the case of the three- 
dimensional wing because more attention has been given to the former 
by other authors. (See references 1, 2,. 7 V 8, 9,. and 10.) It will 
now be shown that the blockage correction is again given by 
equations (19) and (20) where t has a value appropriate to a cir- 
cular tunnel. Since K 3 and K 4 depend on the model and not on 
the tunnel, they are still given by equations (22) and (23). 

The most convenient starting point is the for mula of 
references 1 and 2, namely. 


Ml „ tx/SbV '' 2 
u- \cj 


(25) 


where Sm is the maximum cross-sectional area of the model, 
only necessary to note that 


It is 
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S m = nt 2 /k 
V = K 2 ct2 

It follows at once that 


AiU' , n 3 / 2 t 3 \t n 312 ct2 

— = t_ 8 "c 375 “ t T ” 5 ” c 3 ^ 

/ 


(26) 


Substitution from equations (22) and (23) reduces this to 
equations (19) and (20), respectively. It should he noted that t 
of references 1 and 2 is identical with t . of the present report. 

The value of t is given in table I and figure 2. Values of 
K 3 and K 4 are given in tables IV and V. Again equation (19) is 
preferable to equation (20) . . 

The result of reference 7 fails to take into consideration the 
body shape and so it is not very useful except, for less exact cal- 
culations. The result of references 8, 9, and 10 is presented in 
a different form, namely,, 


where 


AiU' 

U* ' 


1 T ' V 

X V T V 53 


\y factor depending on model shape 

Ty factor depending on tunnel shape 

V volume of model 

B diameter of wind tunnel 


(27) 


It is of interest to show that equations (19) and (20) can be deduced 
also from equation (27.) showing the latter ,to be equivalent to 
equation (25). From reference 8, there is obtained 
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Slime also C = jtB 2 /4 equation (27) yields at once 

AiP' _/« t 3 , N/lt „\ ! V . U eta 

u» ^7 c 3 / 2 y c '8 cV2 

This is the same as equation (26) which yields equations (19) and 
(20) directly. 


Three-dimensional wing .— The blockage correction for a three- 
dimensional- wing in a circular tunnel is given in references 8 , 9, 
and 10, the formula being of the same form as for a body of revo~ 
lution in a circular tunnel, namely, equation (27) where ry and 
\y fcave values appropriate to the three-dimensional wing. From 
reference 8 (using the notation of reference 3 instead of refer- 
ence 2) there is obtained 




V 


= it T 


Since also C = JtB 2 /4 equation (27) yields at once 


(28) 

(29) 


AxU' 

U» 


g sc 2 A Vit tVA 372 1 \_ 2sct A !t3/2 

8 V A*' J \ 8 C 3 / 2 J~ C 3/2 t/c 16 


This is the same as equation (7) which leads directly to 
equations (8) and (9) with the same definitions of K x and: Ks 
given in equations (ll) and (12); T is, however, given by 
equation (29) where Ty is obtained from references 8, 9, or 10. 

Thus equations (8) and (9) may also be used for circular 
tunnels provided only that the appropriate values of r are used. > 
Values of t are given in table I and figure 2. Values of Kx 
and Ka are given in tables II and III. 
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SOLID BLOCKAGE IN .COMPRESSIBLE FLOW 

In the preceding section the soli d-h lockage corrections have > 
Been determined under the assumption that the fluid is incompressi- 
ble. It is now necessary to determine the modifications required 
in these formulas to take account of the effects of the compressi- 
bility of the fluid. The methods of references 12 and 13 are very 
convenient for this purnose. As the required modifications are 
given incorrectly in references 5 and 6 partially incorrectly in 
reference 4, and correctly in references "7> 8, 9 ; and 10, it 
appears worth while to give some discussion of the matter. 

For the purpose of deducing, the properties of a compressible 
flow from those of a corresponding incompressible flow the so-called 
"Extension of the Prana tl Rule," which was first given in reference 12 
and repeated as Method IV in reference lj-, is probably of most 
general application; but the other methods of reference 13 are 
sometimes more convenient for certain problems. The Extension of 
the Prandtl Rule may be expressed in the following manner: 

The streamline pattern of a compressible flow to be calculated 
can be compared with the streamline pattern of an incompressible 
flow which results from the contraction of the y— and z— axes includ- 
ing the profile contour by the factor (M = free— stream 

Mach number) (x— axis in the direction of the free stream). In the 
compressible flow the pressure coefficient as well as the increase 
in the longitudinal velocity are greater in the rat io l/fl-M^) 
and the streamline slopes greater in the. ratio l/^/i-M 2 ■ than those 
at the corresponding points of- the equivalent incompressible flow. 

Since formulas (8) and (19) for the three-dimensional wing and 
the body of revolution have the same form and also apply to both 
rectangular and circular tunnels, it suffices to determine the 
modification due to compressibility for them. Let the subscript c 
refer to compressible flow and the subscript i .to the correspond- 
ing incompressible flow. If V c is the volume of the model in the 
compressible flow, then V±, the volume of the model in the 
corresponding incompressible flow, is given by 

Vj = [l-(M') 2 ] V c 


where M' is the apparent free -stream Mach number at the model as 
determined from measurements taken at a point far ahead of the 
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model. Also if Cc is the cross-sectional area of the tunnel in 
the compressible flow, then the cross-sectional area of the 

tunnel in the incompressible flow is given by 


Ci = il-(M') 2 ] C c 


t is unaffected by the transformation and the effect on Ki and 
K 3 is sufficiently small that it may be neglected. From the 
Extension of the Prandtl Rule it follows that 



1 ( AjU 1 ^ 1 K.jTVj _ 1 K.jTVc 

1-(M' ) 2 U» \ ~ 1-(M’ ) 2 Ci 3/2 [1~(M> ) 2 ] 3 / 2 C c 3/2 


where j denotes the numbers 1 or 3 . .Since equations (8) and ( 9 ) 
are equivalent as are also equations ( 19 ) and ( 20 ), it follows that 
in all cases it is only necessary to multiply the blockage correc- 
tions given by these formulas by ) 2 ] “V 2 in order to take 

account of the compressibility of the fluid. 


As a check it is useful to determine the compressibility mod- 
ification for the case of a body of revolution in a. circular t unn el 
by Method II of reference 13, since the derivation is so’ simple by 
this method. Both the body shape and the longitudinal velocities 
are the same in the corresponding compressible and incompressible 
flows; onl y the tunnel dimensions are altered by the factor 
a/I— (M’ ) 2 so that Ci = [l— (M 1 ) 2 ]Cc- There is obtained 



_ &a T Vl _ 1 K 3 ty c 

C^ [l-(M' )2 ]3/s c 3/2 

c 


as before. 


WAKE BLOCKAGE IN COMPRESSIBLE FLOW 

The blockage due to the wake of a model in a two-dimensional- 
flow tunnel is discussed in some detail in reference 3 . Much of 
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this discussion is applicable without change to the case of a three- 
dimensional— flow tunnel. The fundamental idea of replacing the 
model and its wake by a source (in this case a three-dimensional 
source) of suitable strength located at the position of the model 
can be used again and in fact the determination of the source 
strength 0 can be carried out in exactly the same manner. The 
result is identical with that of reference 3* namely. 


Q = ^D -. S . |1 + ( 7 _1)(M')2] (30) 


where 


Q. mass flow of source rather than volume flow as used previously 

p* mass density of fluid at point far upstream 

•V • uncorrected drag coefficient referred to apparent dynamic 
pressure q* 

S area on which drag coefficient is based 
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ratio of specific heat of 
heat at constant volume 



constant pressure to specific 


In equation (30) powers of M* higher than (M* ) 2 have been 
neglected. 

Consider now a rectangular tunnel of height H and breadth B. 
The tunnel walls are replaced by 3. doubly infinite array of. sources 
of strength Q at distances mB to the side and nE above and 
below the position of the model.. 


By making use of method II. of reference 13, it is readily shown 
that a three-dimensional source of strength Q. (mass per unit time) 
in a uniform flow of compressible fluid will induce at the point the 
coordinates of which are x, y, z relative to the source a stream- 
wise velocity 


Au 


_ _JJ X 

^Jto [x 2 + (l- M 2 ) (y^z 2 )] 3 / 2 
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where the uniform flow is in the x— direction and p and M are the 
density and Mach number, respectively, of the undisturbed stream. 

It follows that the streamwise velocity A2U 7 induced at a 
point on the center line of the tunnel by the entire system of 
images is 


AaU« = JL 

^itp* 



x 

[1-(M‘ )2 ] (m^+r^H 2 ) | 3 / 2 


where p ’ and M* are the density and Mach number of the undis- 
turbed flow in the tunnel and the summation is taken over all 
positive and negative integral values of m and n except 
m = n = 0. The velocity induced by the image sources at an infi- 
nite distance upstream is ; 


(A 2 U f ) = lim -Q — \ 

t co x oo 4np' ^ x 2 + [1— (m* ) 2 ] (m%s 2 +n 2 H 2 )'l 3/ '’ 2 


But conditions far upstream must remain -unchanged and to achieve 
this it is necessary. to counterbalance this velocity by super- 
position of a uniform flow of equal magnitude but opposite sign. 

The addition of this flow at all points in the field will result in 
a speeding up of the general flow at the position of the airfoil by 
the amount • 


A 3 U’ 


lim 

oo 4rtp 


X 

jx 2 + [1— (M* ) 2 ] (m 2 B 2 +n 2 H 2 ) \ 3 / 2 


(31) 


The summation of this series can be obtained by the following 
artifice. .Substitution of equation (30) into equation (3I) and 
setting M’ = 0 gives 



= lim 

X -> 03 


Cp'SU 1 
8 n 


f 

\ £ ; 

/ s ( x 2 +m?B 2 +n 2 H 2 ) 3 /2 


( 32 ) 
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But according to references 4, 8 , 9 , and 10 


Agr \ _ 1 Cp'S 

. U* J t ~ 4 BH 


(33) 


Comparison of equations .(32) .and ( 33 ) shows that 


lim V — = !=* ( 34 ) 

x 00 [_J ( x 2 +m 2 B 2 +n 2 H 2 ) 3 / 2 BH 


If in equation (34) B and H are replaced By )2 B 

) 2 respectively, there is obtained 


and 


lim y 

X ■» OB / . 


• X 

j" x 2 +£l— (M* ) 2 ](m 2 B 2 +n 2 E 2 )\ 3 / 2 


2rt 

[1_(M»)2]BH 


(35) 


Substitution of equations ( 35 ) and ( 30 ) into equation ( 31 ) yields 


AgJ 1 _ 1+(7-1)(M ') 2 Cp'S _ 1+0. 4(M* ) 2 Cp’S 
V 1-(M' ) 2 4BH l-(M ') 2 4C 


on setting 7 = 1.4. The preceding discussion is for a rectangular 
tunnel, but in references 8 , 9 , and 10 the same formula is given for 
a tunnel of any shape for the case M* - 0. Consequently, equation ( 36 ) 
may be taken to hold generally. for a tunnel of any shape. 

In reference 4 the wake— blockage correction is given correctly 
for incompressible flow but the modification to take account of 
compressibility is given incorrectly, as it is determined from ref- 
erence 5 which," as pointed out in reference 13 , is incorrect. The 
effect of compressibility on the wake— blockage correction is deter- 
mined in references 8 , 9 , and 10 by means of the correct form of the 
Extension of the Prandtl Buler^but there is some question whether this 
rule is applicable to this problem. It appears that in using this 
method to determine the effect of compressibility on the blockage 
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correction, the effect of compressibility on the source strength Q 
given by equation ( 30 ) is overlooked. This explains the discrepancy 
between equation ( 36 ) and the.fomula of references 8 , 9 , and 10 .' 

WAKE PRESSURE GRADIENT IN COMPRESSIBLE FLOW 

B 'P 1 ? effect of the pressure gradient caused by the wake is dis- 

cussed in reference 3 for compressible flow in a two-dimensional 

wind tunnel and in reference 11 for compressible flow in a circular 

the The + Qethod of reference 11 is equally applicable to 

T ^ a - reCtangUlar vind tunnel if the appropriate value of 

T stable I; is used. 

+ nrm P 1 ® longitudinal velocity increment due to the effect of the 

a3S-i^^ ar 1 ie3 ° n the SOUrCe USed to slmulat3 the wake has an 
loca^iT sradient in the stre am direction at the model 

lineer ^?^ gradient in the velocity, is equivalent to a 

mte^eStlon “ “ S 18 eaBliy 300n fr0E tho ^PProxi- 

Ap' — 2q< ^21 

In reference 11 it is shown that the gradient in AsU* for a 

f nd tUme i f B ident ically equal to the value of AeU' 
for a doublet m a wind tunnel. In the notation of the presen? . 
report there is obtained P 

(&-) =- o' JTt C D' S 

'dx ' . ^ h c 3/2 


For compressible flow this becomes 


dp 1 

dx 


= — q* ±-± (7- 1)/M') g . J~n Cp'S 


~ -9 


* 1 + O.h 


[ 1 - (M» 




. <in Th ® increase in the drag resulting from this pressure gradient 
equal to the product of the pressure gradient by the sum of the 
actual model volume and the virtual volume (reference 2) It 
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should be noted that the constants Kj and Ks of tables II and IV 
are equal to ~j\l ) .times. the ratio of the sum of the actual model 

volume and the virtual volume to the actual model volume, these 
quantities being calculated under*~the assumption of incompressible 
flow. Thus the increase in drag coefficient caused by the pressure 
gradient is given by 


ACd w ’ 


i + o.w) 2 K i*yy 

[1 _( m «) 2] 3 /2 C 3/2 


(37) , 


for the wing . and 


ACpb' 


1 + Q.kfM 1 ) 2 K 3' r yb C D * 
[1 -(M') 2 1 3/2 c 3/2 


(38) 


for the body of revolution. 

It is pointed out in reference 11 that the virtual volume is 
altered by the compressibility of the fluid. Thus equations (37) 
and (38) can be slightly improved if Ki and .Kg appearing therein 
are corrected to take account of this effect. The '-necessary modifi- 
cation is made by replacing Ki and K 3 in these equations by 
Kip and Kbp, respectively, where 



5 = 1,3 

' 


(39) 


where h(M' ) is given by figure h which is reproduced from 
reference 11. It should be npted that the . improvement in modifying 
Ki and K 3 for compressibility will be ' small especially in the 
case of K 3 for the body of revolution. It is important only for 
quite high Mach numbers. , This additional ref inement was not made 
in reference 3* 

It appears that a similar compressibility modification of Ki 
and Kq in l 116 formulas for the solid blockage should be made, but 
the exact modification required is not known. However, it is 
believed to be small. 
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CORRECTION OF MEASURED QUANTITIES 


The true velocity U at a model consisting ^ 'body of revolu. 
tion and wing can he obtained from the apparent velocity u * 
ap^lyiSg the solid-blockage and wake-blockage corrections. The true 

velocity may be written in the form 


U = U’ (1 + K) 


CM» 


where 


K = Ky. + 1^, + 


(4l) 


In equation ( 38 ) K„ is the solid-blockage correction due to the 
wing and is given by 


or 



[ i -( m ») 2 ] 3/2 


Ki t Vw 

c 3 /2 


1 • KgT2sc v .t v 

= [1-(m') 2 1 3 ' s c 3 ' 2 


(42) 


(43) 


is the solid-blockage correction due to the body of revolution, 
being given by 0 


! V^b 

Kb “• [ MM - c 3 / 2 


(WO 


or 

i 1 ^4 Tc b^b 2 (45) 

K b " [ 1 — ( M * ) 2 ] 3 / 2 C 3 / 2 

is the ' wake-blockage correction and is given by 
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. 1 + O.^M 1 ) 2 Cp'S 

1— (m* ) 2 4c 


It is evident that a correction to the apparent velocity in a 
compressible flow implies corrections also to the apparent density, 
dynamic pressure, Reynolds number, and Mach number. These correc- 
tions are readily obtained on the basis of the usual assumption that 
the flow is adiabatic. It is assumed that the correction terms are 
small compared with unity, so that squares and products of these 
terms may be neglected. The analysis follows the lines of reference 3 
and 'it is not necessary to repeat the details here. The following 
equations are obtained: 



p = p-*[i - (m») 2 k] 

(^7) 

q = q’| 1 + [ 2 — (M* ) 2 ]K 1 

(W) 

R =R'-|l + [l— 0.7(M* 

(^9) 

1 

M = M’-f 1 + [1+0. 2 (M* ) 2 ]kT 
l ■ J 

(50) 

/ 

The drag coefficient must be corrected for the effect 
pressure gradient due to the wake as well, as to refer 
correct dynamic pressure. There is thus obtained 

of the 
it to the' 

C D = V j 1 ~ C 2 -(M’ ) 2 3K - ACpJ - AC^ } 

S 

(51) 


where ACp^ and ACd^* are given by equations ( 37 ) and ( 38 ). - 

Numerical values of the functions of M* which appear in these 
equations are given in table VI. 


CONCLUDING REMARKS 

Data obtained from tests of three-dimensional models, which are 
small relative to the wind-tunnel dimensions, can be corrected for 
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solid and wake 
wake by means 

blockage and' for the pressure gradient due to 
of the following equations: 

the 


U = U* (l+K) 

(ito) 


q = q' |l + [2-(M’) 2 ]k| 

(48) 

i 

R = R' jl + [1-0.7(M’) 2 ]k| 

(49) 


M = M* “jl + [1+0.2(M') 2 3 k| 

(50) 

C D 

= C D ’| 1 - [2 ~(M« ) 2 ] K - AC Dw » - j . " 

(51) 

In the preceding equations K is obtained from: the following: 


K = K w + K b +K wk 

(41) 

where 

1 K i TY w 

. [1-(M’) 2 ] 3/2 e- 3 / 2 

(42) 

or 

y ■_ 1 K 2 T2sc v t w 

■ " [1— (M* ) 2 ]3/2 C 3 / 2 

(43)' 


r . _ 1 Wb ' • 

° [l-(M* ) 2 ] 3 / 2 C 3 / 2 

(44) 

or 



/ 

Y - 1 K 4 1 ' c b t b 2 

* ~ [ 1— (M* ) 2 ] 3 / 2 C3/2 

(45) 


. _ 1+0.4(M') 2 CD' 5 
wk 1-(M* )2 4C 

(46) 



MCA EM No. A7B28 


24 


AC Dw' 


1 ± 0. 4(M* ) 2 KjT W 
[1 -(M*)^ 3 / 2 C ®/a 


(37) 


AC , - 1 1 0.4(MM S KstVdCd 1 

^ [1-(M') 2 ]3/ 2 .C3/2 


(38) 


In these equations t is a factor which for a body of revolu- 
tion depends only on the shape' of the tunnel; whereas for a three- 
dimensional wing t depends on the ratio 'of ' the wing span to the 
tunnel breadth as well as on the tunnel shape. (The wing span is 
in the direction of the tunnel breadth.) Values of t are given in 
table I and figure 2. The constants Ki and Ks for the three- 
dimensional wing depend only on the wing base profile shape and can 
be calculated by means of equations (11), (12), and (13)- Values of 
Ki and Kz for a number of wing profiles are given in' tables II and 
III. The constants K 3 and HU for the body of revolution depend 
only on the shape of the body and can be calculated by means of 
equations (22), (23), and (24). Values of these constants for a 
number of body shapes are given in tables IV and V; some of these 
shapes are drawn in figure 3* Since these tables are incomplete and 
fuselage forms are not standardized as are wing sections, it is 
recommended that the values of K3 given in table IV for the 
NACA 111 series of shapes be used for any fuselage shape which does 
not differ too greatly from an NACA 111 shape. This implies that 
equation' (19) and table IV should be used in preference to 
equation (20) and table V whenever possible. Numerical values of 
the functions of M' which appear in the correction equations are 
given in table VI. 

The constants Ki and K 3 appearing in equations (37) and (38) 
may be modified for compressibility by means of figure 4. The 
modified values of Ki and K3 are to be used in equations (37) 
and (38) only and not in equations (42) and (44). 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, • 

Moffett Field, Calif. ■ . 
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APPENDIX A 

SUMMATION OF THE INFINITE SERIES FOR a (0, s/B, B/H) 

From equations ( 3 ) and (16) it is seen that v (0, s/B, B/H) 
may he written in the alternative forms 


( 0 } §. b \ iV'l ( m±gjB _ ie-s/B ( A1 ) 

V ’B'Hy 2s Lx0 V V F^ J 


or 


N V 


q s B . \ _ \ • 2m • 

^ ^ y *~ J EumFinnt (m+ s /B ) F mn+ (m— s /B .) Enin] 


'(AS) 


where the summation is taken for all positive and negative integral 
values of m and n except m = n = 0 and the quantities E^ 
and Fjgjp which are introduced for convenience,, are defined by the 
equations 


Emn = v n 2 +(m+s/B) 2 (B/H) 2 
■^mn = S n2 + (m — 3 /B ) 2 (B /H ) 2 


It is possible to sum the series for 0 (0, s/B,' B/H) exactly 
only when s/B = 3 . In this case equations (Al) and (A2) yield 



26 


NACA m No. A7B28 




FP 

eg 

I — I 

I 


8 P-PTh] 

OJ 

+ 

I 1 


H| 0 j 

I 

a 


Hjoi 

+ 

a 


S 

I 

s 

&k 

w 

p=r 

eg 

hTw 


i 


eg 

fi]8? 

I 


1 


OJ 


rHiCVJ 

IXlU 


OJ 


OJ 


PQ 


eg 


FP 

eg 

rjM 

vE 

.+ 

w g 

\ 


8 IP J a 

a 8 

•H A 

a . 

H ^ 

leg 

|0J 

h| a 

h la 


8 

a 

OJ 

OJ 


> i — 1 

1 h 


i 1 

, eg 

j? 
FP 


«|VD 

I J. 

FP 

CO^ 
FP 


h|>& 

•Ki 

FP 


H 

II 

a- 


CO 
■ © 

•H 

Jh 

© 

CO 

© 


H I 
gJ 
o 

•H ' 
U 

1 Q 


© 

© 

•H 

JP 

© 

CD 

© 

A 

-P 


O 

*P 

u 

e3 

© 

© 

© 

O 

© 

a 

© 

•ri 

•P 

*H 


JW 

co 


Ch 

O 

© 

• © 

PS 

r H 

03 

> a 

© 

a +> 

© -P 

rCf *H 

S .Si 

© 

u u 

o 

Fr o 

fQ 


eg 

© 

+ 

a 

s. 

w « 


OJ 

+ 


w 

5 

GO’ 

A 

+ 

Fr 

PP 

© 


a 




G 



on 



s u 


OJ 


(B/H) 3 / tm 2 -(s/B) 2 ] 



where N is an arbitrary positive integer and B N (0, s/B, B/H) denotes a remainder term, 
evaluating a(0, s/B, B/H) in the present report, N was usually taken equal to 7 and the 
summations indicated in equation (A 3 ) were carried out exactly; whereas an approximate value 
was used for the remainder term. An approximate formula for Bn( 0, s/B, B/H) is 
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provided a, b >0. Now K N (0, b/B, B/H) is approximately equal to 
0, B/H) and it is easily found that 

00 00 00 N 

%(°' 0, B / H) = 4 ( Y, X + X X ) [ n 2 +m 2 (B /H ) 2 ] 3 / 2 

n=l m=N+l n=N+l m=l 


00 ' 

. ? V i + _J 

L n 3 (B/E) 3 

n=N+l 


) 

L j 

m=N+l 


m c 


These summations may be approximated by suitable integrals so that 
approximately 


% 


r 


n s B 

i '-'j — > 

\ B H 


=k 


00 


. n" 

dy / 

JN+ A 



r>°° 




1 

[ y 2 + (B /h) 2 x a ] 3 / 2 


+ 2 




(A7) 


If the integrals of equation (A7) are evaluated by means of 
equations (A5) and (A6) , then equation (At) is obtained. 
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APPENDIX. B "■ 

LIST OF IMPORTANT SYMBOLS 


B tunnel "breadth or diameter 

H tunnel height 

C tunnel cross-sectional area 

c chord of airfoil or body of revolution 

.t maximum thickness of airfoil or body of revolution 

V volume of wing or body of revolution 

S m maximum cross-sectional area of body of revolution 

s half span of wing 

Cp drag coefficient 

S model area on which drag coefficient is baaed 

U stream velocity • 

M Mach number 

R Reynolds number 

.7 ratio of specific heat of gas a.t constant pressure to 

specific beat, at constant volume (cp/c v ) 

p mass density 

q dynamic pressure 

A,KijK 2 factors depending on shape .of airfoil base profile 
(See equations (ll), (12), and ( 13 ) and tables II 

and III* ) 

factors depending on shape of body of revolution (See 
equations (22), ( 23 ) and (24) and tables IV and V.) 

linear compressibility correction factor for virtual 
volume (reference 11 ) 


h (M * ) 
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factor depending, on tunnel shape and wing— spa i>- to— 
tunnel— breadth ratio (See equations (10) and (21) 
and table I.) 

total blockage correction (See equation (4l).) 
wing— blockage correction (See equations .(42) and (1+3) • ) 
body— blockage correction (See equations (44)and (45).) 
wake— blockage correction (See equation (46) . ) 

Superscript 

( 1 ) when pertaining to fluid properties, denotes values 

existing in tunnel far upstream from model; when 
pertaining to airfoil characteristics, denotes 
values in tunnel, coefficients being referred to 
apparent dynamic pressure q' ' 

Subscripts (used only when necessary to avoid ambiguity) 
c denotes values in compressible fluid 

i denotes values in incompressible fluid 

w denotes values for wing 

b denotes values for body of revolution 
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K 

Kw 

K b 

Kwk 


wk denotes values for wake. 

i 
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TABLE III.- VALUES OF K 2 FOR VARIOUS BASE PROFILES 
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TABLE VI— COMPRESSIBILITY .FACTORS FOR CORRECTION EQUATIONS 
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Values of t for various tunnel shapes and conf igu rations. 
















NACA RM No. A7B28 


Fig. 2b 
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SPAN-TO-TUNNEL- BREADTH RATIOS. 
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Fig. 4 
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